The Topological Properties of Fault-Tolerant Circuits 



O 

o 

(N 

> 
O 

(N 



> 

0^ ■ 
O ■ 

(N ■ 
O 

o ■ 



Alastair Kay 

Centre for Quantum Computation, DAMTP, Centre for Mathematical Sciences, 
University of Cambridge, Wilberforce Road, Cambridge CB3 OWA, UK 

The design of time-independent local Hamiltonians that realise quantum algorithms is derived 
from the study of perfect state transfer. The novel features of this evolution are the perfect realisation 
of the computation, and the ability to implement operations in parallel rather than sequentially. 
This is applied to a quantum circuit for concatenated error correcting codes, and the properties 
of the eigenstates are compared to those of topological memories. An erroneous implementation 
of the algorithm (one where a single error could lead to a fault in the output) is found not to be 
topological, whereas a true fault-tolerant algorithm possesses self-correcting properties, suggesting 
the existence of a unifying structure for fault-tolerant memories. 



Introduction: During a quantum computation, it will 
be vital to protect the system from noise, and to min- 
imise the effect of any residual decoherence with the help 
of some form of error correction or decoherence-free sub- 
space. The first important step towards this is a quantum 
memory, which would provide us with the ability to store 
a quantum state, protecting it against decoherence, for 
long times. There are two very different ideas for imple- 
menting such a memory. The first, concatenated fault 
tolerance (CFT), uses a series of concatenated error cor- 
recting codes m 0, 13 • The second idea is to implement 
a topological quantum memory (TQM), where quantum 
information is stored in the degenerate ground states of a 
local Hamiltonian. To enable storage of information for 
an arbitrarily long time, schemes such as the toric code 
Q are not sufficient as they do not protect against ther- 
mal noise. Instead, self-correcting codes are required (a 
system is self-correcting if the survival time of the stored 
information ~ e"^/'^ with system size N , when subject 
to local, thermal noise of temperature T), of which there 
is a known example in four spatial dimensions 0], and 
one is conjectured in three dimensions [5[. 

In this paper, we make steps towards a unifying struc- 
ture for these two types of quantum memory by con- 
structing a fixed Hamiltonian whose dynamics implement 
a CFT scheme, and is spatially local on a 2D lattice. The 
eigenstates of the Hamiltonian correspond to a superpo- 
sition of all the steps through the computation, and is 
degenerate, with the degenerate space being described 
by a set of orthogonal input states. By adding a series 
of energy penalties to the Hamiltonian, the degeneracy 
can be lifted such that only the logical states of the error 
correcting code are in the degenerate ground state space. 
Once in this position, we are able to discuss the topolog- 
ical properties - how easy it is for local noise to perform 
an SU (2) rotation on the ground state space. 

The simplest approach to examining the similarities 
between CFT and TQM would involve taking the defini- 
tion of an error correcting code in terms of stabilizers, and 
then write a Hamiltonian as a sum of stabilizers. This 
approach doesn't interest us, primarily because the resul- 
tant Hamiltonian is not local, and as soon as we use non- 



local terms, we ought to consider the effect of non-local 
errors etc., and we lose all comparison to TQM. A second 
consideration is that our approach allows us to discuss 
the implementation of a fault-tolerant circuit. This is a 
very important aspect since we know that the construc- 
tion of schemes that do not propagate errors catastroph- 
ically is both necessary and difficult. One might expect 
a phase transition between the cases of perfect (pCFT) 
and imperfect (iCFT) implementation of CFT and, in- 
deed, we observe this feature. 

Our derivation of the Hamiltonian follows on from the 
study of perfect state transfer, a popular subject of re- 
cent years (see, for example [y, 0, S 0), which shows 
how a fixed Hamiltonian can be designed that perfectly 
transfers a single excitation from one end of an exchange- 
coupled chain of qubits to the other. These ideas have 
already been applied to other situations where we wish 
to transfer from one state to another, such as the im- 
plementation of a controUed-NOT gate and quan- 
tum state amplification 11]. The ultimate generalisa- 
tion of this, which we demonstrate here, is the design of 
a Hamiltonian to evolve from the input of a computa- 
tion to the output. The construction is closely related to 
the ideas of Fcynman [l5| , and Kitaev (used to show 
that the problem of finding the eigenvalues of Hamiltoni- 
ans is QMA-complete). In comparison, our Hamiltonian 
achieves the desired computation perfectly, and incor- 
porates a degree of parallelism. Perfect state transfer 
schemes tell us that if there is a set of orthogonal states 
{| "(An)} which are acted on by a set of Hermitian opera- 
tors {i?„} such that 

Rn I V"™) = 5n,m \ V'n-l-l) + ^n+l,m \ V'n) , 

then there exist sets of coupling coefficients {Jn} such 
that the Hamiltonian H — 'Y2n=o JnRn implements the 
transfer | i/io) — > | V'Af) in a fixed time, t^. One example 
is given by J„ = \/ (n -|- 1)(A^ — n), although there are 
many other analytic Q and numeric solutions [3| . 

Quantum Computation: Let us consider a general com- 
putation of N distinct steps J7„, neglecting the measure- 
ment step, which we shall assume occurs at the end of 
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the computation, since if a part of the computation re- 
quires the result of a measurement, it can be performed 
as a controlled gate instead. The state after n steps of the 
computation, | = t/„ . . . J7i | iPq), is not necessarily or- 
thogonal to the states from the previous steps. This can 
be avoided via the introduction of a clock, whose value 
starts at |0)p, and is incremented by 1 as each gate is 
applied, so | "ifn) ^ \ n) ^ \ tpn) , (*m|^'ri) = Snm and 

Rn = \n + 1) {n Ip (g) Un+i + \n) {n + l\^(g) ul^^. 

To implement the clock, we use a set of -I- 1 qubits, and 
the states | n)^ denote the location of a single excitation 
within the set, giving 

Rn = (o-,;ct++i)c «) Un+i + (cr,t<+i)c «) J/I+i 

and the Hamiltonian is implemented with only 4-body 
terms. The computation is initialised with the computa- 
tional qubits in the input state, and the clock in the | 0)^ 
state, and the finishes when the excitation is detected on 
clock qubit TV, which is guaranteed at time to. Following 
14j . this construction could be reduced to fewer terms. 

Implicit Clocks: If we wish to demonstrate fault- 
tolerant properties in a Hamiltonian, the whole system 
will need to be fault-tolerant. However, if we use a 
scheme involving a clock, we would somehow have to 
build in protection for the clock as well as the circuit. 
Instead, we shall implicitly involve a clock by using spa- 
tial position, in a similar fashion to the ID scheme of 
Gottesman et al. [3]. Let us consider a square lattice 
in 2D (with periodic boundary conditions), composed of 
qutrits. The whole system will be initialised in the state 
I 2) , except for a single row, which is the input state of 
the computation, defined on the | 0) and 1 1) levels of the 
qutrits. The clock of the computation will be implicitly 
defined by the location of the 'active' region of the com- 
putation - this will proceed by moving the | 0) and | l)s 
from row to row, one qubit at a time. For example, if 
on row 2n, cell m, the qutrit to the left is either | 0) or 
1 1) (i.e. not I 2)) and the qutrit to the right is | 2), then 
we should swap the qutrit at that position with the one 
on row 2n — 1. On rows 2n -|- 1, we swap left and right. 
In addition, we should implement a unitary operation J7' 
between the qutrit being moved and one of its neighbors. 
Overall, there is a local 5-body Hamiltonian term 



= (l-|2)(2|)2„.™-i |2)(2|,„^„_,, 5g:i™!™)C^S" 



l,m) 



where is a swap operation between qutrits i and j. 
Terms at the ends of rows are adapted such that the 
evolution continues to propagate (i.e. the active region 
should be propagated forward two rows). 

Parallelism: So far, our construction only involves a 
single localised region where computation is proceeding, 
whereas Aharonov et al. have shown that a parallelism 
of at least log(iV) is required for fault-tolerance [3|]. The 
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FIG. 1: A generalised hypercube construction. The depicted 
case is two-dimensional, with A'^ + 1 qubits along each direc- 
tion. This can be treated as an effective spin chain, as shown. 
When used in computation, the two dimensions of the hyper- 
cube correspond to two parallel threads of computation. 

study of perfect state transfer contains sufficient insight 
to incorporate the desired parallelism. For state transfer 
between antipodal points on a fc-dimensional hypercube 
Q, the excitation travels along all fc-dimensions inde- 
pendently, precisely as required. There are two simple 
additions that are necessary. Firstly, each of the parallel 
threads should be the same length. This is easily realised 
by padding the shorter threads with 1 operations. The 
second is that we need to incorporate the hypercube into 
our larger state transfer structure. This is achieved using 
the constructions of [10]. Specifically, state transfer on a 
hypercube can be written as an effective transfer along a 
spin chain, as depicted in Fig. [TJ If it is organised such 
that the effective chain has those couplings that we are 
using the hypercube to replace, then the system cannot 
tell the difference, and state transfer occurs as before. 

Multiple working areas can be prepared by using the 
same mechanism that we use to reverse the flow of a com- 
putation - to shift some cells further up than you other- 
wise would. Provided that parallel threads act on spa- 
tially distinct areas of the computer, the computer func- 
tions as desired, and when we wish to recombine parallel 
threads, they all do so before the rest of the computation 
proceeds, by virtue of the fact that the active region of 
the computer is unable to move until all computational 
qubits in the neighborhood are present, i.e. all the local 
threads have terminated. 

It is beyond the scope of this work to give explicit ex- 
pressions for coupling strengths on the hypercube which 
achieve the desired effective spin chain. In principle, 
there is sufficient freedom available - for a cube of dimen- 
sion fc and A^-|-l qubits along each direction, there are kN 
couplings Jn and also fciV couplings in the effective chain. 
Unfortunately, the equations incorporate terms such as 
I J„P, which must be positive, and hence there is no guar- 
antee that solutions exist. In general, numerical solu- 
tion is required to find useful coupling schemes. In order 
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to justify that solutions exist, we will demonstrate one 
scheme that is of particular relevance, where the effective 
chain is uniformly coupled, in the special case of k — 2. 
We shall take the couplings along one dimension to be J„, 
and along the other to be Kn = Jn-i-u- The symmetry 
guarantees that the effective coupling strengths satisfy 
L2n~i = Ln, which are given by 



Jr + Kr 



(1) 



where Ji — YVj=o 1"^!^ ^^'^ < r < 2N. Perhaps the 
most illuminating term is the central one, where 



L 



N 



~ 2^+\2N+iy.l 



A simple way of satisfying this equation is to set Ji ~ 
2(^^ij 1 which solves all Eqns. ([T|) to a good approxima- 
tion; while the first term Lq = 3/4, the second term, in 
the limit of large N, is = 35/36, and the accuracy 
improves towards the centre of the chain. This accuracy 
can be further improved by slight modifications of the 
strengths at the ends of the chain. Note that the con- 
struction for higher dimensions can be performed in an 
inductive fashion - we create a two-dimensional hyper- 
cube with effective coupling strengths J„, and couple it 
to another single chain with strengths Kn- As a result, 
we have a three-dimensional hypercube, which has an ef- 
fective chain which is uniformly coupled. Furthermore, 
all of the required coupling schemes are approximately 
constant, so one expects that small modifications of the 
uniform scheme should yield good approximations. 

CFT and TQM: We are now in a position to implement 
a CFT circuit, and explore how the resultant Hamilto- 
nian relates to topological models for quantum computa- 
tion. We shall take a 2D x M lattice. Across the width 
of the lattice, we encode N/n logical qubits (in the first 
level of concatenation) of distance 3 [l8[, which will im- 
plement R levels of concatenation, so = n^. The depth 
of the lattice M ~ poly(Af) corresponds to the number 
of time steps in the circuit. Fault-tolerant schemes in- 
volving nearest-neighbor interactions are already known 
liL 17 1 . The model that is closest to that required for 
the Hamiltonian construction is from [l^, which incor- 
porates nearest-neighbor interactions, coherent feedback 
from ancillas and even the concept of read-heads, local- 
ising the active region of the computation (although the 
precise layout of the computer is slightly different) . We 
shall assume that we are using such a scheme, and that 
it has a fault-tolerant threshold of ec per qubit per gate. 

In order to be able to compare our Hamiltonian to 
topological Hamiltonians, there are two technical details 
that must be addressed. The first is that, for our model to 
be reasonable, we must provide protection against errors 
on the entire system, not just the subspace spanned by 
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FIG. 2: Error correction of the 3-qubit bit-flip code, using a 
coherent feedback process. The construction is not suitable 
for concatenation for a CFT scheme since the indicated X 
error propagates to a second position, so after another round 
of error correction a logical bit-fiip has occurred as a result 
of a single error. The circuit can be altered to deal with this 
effect, and it is these two types of circuit that we contrast 
within the framework of Hamiltonian evolution, but for the 
storage of quantum information. 

the I 0) and 1 1) states. This can be achieved by applying 
Ising terms |22)(22| on all nearest-neighbor pairs. This 
makes the state where all spins arc in | 2) a self-correcting 
ground state of the system. Of course, this has the con- 
sequence that none of the eigenstates with overlap on the 
initial state of our system are ground states. However, 
it is not an essential property of topological memories 
that they be ground states (although it is preferable). 
Let us define by | Ol) and 1 1^,) the eigenstates of mini- 
mum eigenvalue that have overlap with the initial state 
our system, with the input to the computation being the 
correctly encoded logical | 0) and 1 1) states. These eigen- 
states will take the form 



I Or,) = 



(2) 



i=0 j = l 



where we have suppressed the position information for 
clarity, and there are a total of D steps 19] (for sequential 
operation D — NM). 

The second problem that we must tackle is less satis- 
factory, but is an essential element. In a CFT scheme, 
many ancillas are used for syndrome extraction and cor- 
rection. After use, these are reset, but we are unable to 
reproduce this with Hamiltonian evolution. Since we are 
only performing a single cycle of concatenated error cor- 
rection within the circuit, we only require a finite supply 
of ancillas, which we can therefore take to be present at 
the outset of the computation. However, resetting is still 
a vital component as the ancillas contain a history of the 
errors that have occurred, and hence were we to test the 
signature of topological protection 



(Ol|1^|1l) -0 



(3) 



for two degenerate (orthogonal) eigenstates |0l) and 
I 1l) acted on by a single-site error V, this would nec- 
essarily always hold because of the information stored 
in the ancillas. We must incorporate the tracing out or 
resetting of ancillas when evaluating the condition ([3]). 
Were we to calculate the overlap in Eqn. the only 
terms that could possibly have non-zero overlap are those 
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evaluated at equal times. Thus, we can trace out differ- 
ent parts of the state at different times. Specifically, at 
any particular time step, we can trace out all ancillas 
that should have been reset by that time in the algorith- 
mic implementation. We use this procedure, applied to 
Eqn. (121), as a revised version of the topological condition. 

Consider the case of iCFT, such as the example given 
in Fig, [21 Even if there is only a single location in which 
an error occurs that can lead to the error correction cy- 
cle altering the qubit in the wrong way (to another logi- 
cal state), then within a concatenated scheme, there is a 
single location where that single error in a logical qubit 
causes a flip in the next level of concatenation, and so on, 
to the level at which the real data is stored. Let's make 
this single error the V in Eqn. | 1^) is no longer an 

eigenstate of the Hamiltonian, so we should decompose it 
in terms of eigenstates. Let's consider the component of 
I 1l) which constitutes the time at which V would have 
to be applied for the whole evolution to go wrong. This 
term is part of an eigenvector of the form in Eqn. ^ 
which is not orthogonal to | Ol) because, at the very least, 
the components in the final time step are identical (apart 
from the ancillas, which are traced out). Thus, for these 
non-ideal circuits, we have that (Ol|T^|1l) > l/D. 

In order for there to be an overlap in the pCFT case, 
we need at least ecND errors to occur, and for topo- 
logical effects, there must be an energy gap between the 
degenerate states of interest and the intermediate states 
populated by the errors, which is not currently present in 
our model. Since errors that convert | 0) ^ | 2) can never 
violate Eqn. consider only errors acting in the | 0) and 
I 1) space of the qutrits. Thus, we have restricted to a 
subspace where the number of | 2)s is preserved, within 
which I Ol) and | 1l) are ground states (strengthening the 
comparison to TQM). This is particularly helpful since 
we are able to evaluate the effect of adding energy penal- 
ties to the Hamiltonian thanks to the projection lemma 
of Specifically, by adding terms of strength poly(A^) 
to penalise states that do not conform to our two desired 
input states, corresponding to the logical qubits | 0) and 
I 1) of the CFT scheme, then the relevant eigenvectors 
are increased in energy by at least l/poly(iV). The sim- 
ple way to add these penalties is to utilise the ancillas 
used in the error correction cycles, which detect whether 
the input state was correctly encoded or not. By decod- 
ing the value of the ancillas onto single bits, which can 
be done at a stage where they will no longer affect the 
computation, an energy penalty |1)(1| can be added to 
a single qubit to detect if some correction was required. 
If the energy penalties are all equal, the effect of ther- 
mal noise on the degenerate ground space of | 0^) and 
I li) is similar to the effect on the toric code - as the 
number of errors that occur increases, the energy cost re- 
mains approximately constant (errors build up at 1 level 
of concatenation, up to a fixed maximum, and then re- 
set, propagating a single effective error to the next level 



of concatenation, leaving the total number of effective er- 
rors roughly constant, yielding a constant energy cost). 
Instead, the system can be made self-correcting by vary- 
ing the added energy penalties - we can make the energy 
penalties for ancillas at the m}^ level of concatenation 
to be ^ [d + \y"'J ~ poly(A^) J, so that each error has 
the same energy penalty, J, giving a survival time for the 
stored information of g^cNDj/iksT) ^ 

Conclusions: We have demonstrated how a quantum 
computation can be implemented perfectly using Hamil- 
tonian dynamics. By analogy with state transfer on a hy- 
percube we introduced parallel computation within 
this model. Using this, we have compared circuits for 
CFT and TQM, finding that circuits which are not cor- 
rectly designed for fault tolerance give Hamiltonians with 
no topological protection, because single errors convert 
between the ground states. In contrast, a correctly de- 
signed circuit gives a Hamiltonian with self-correcting 
properties, up to the caveat of a reinterpretation of 
Eqn. ([3l), which was introduced to allow for the reset- 
ting of ancillas within the circuit model. It is our hope 
that this investigation will propel work towards a unified 
structure for fault-tolerance in the differing models. 
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